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The filled-skutterudite compound CeRu4Sbi2 shows a pseudo-gap structure in the optical 
conductivity spectra similar to the Kondo insulators, but metallic behavior below 80 K. The 
resistivity shows a large peak at 80 K, and the Seebeck coefficient is positive and also shows a 
large peak at nearly the same temperature. In order to explain all these features, a simplified 
tight-binding model, which captures the essential features of the band calculation, is proposed. 
Using this model and introducing the correlation effect within the framework of the dynamical 
mean field approximation and the iterative perturbation theory, the temperature dependences 
of the optical conductivity, resistivity and the Seebeck coefficient are calculated, which can 
explain the experiments. 
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The filled-skutterudite compounds RT4X12 (R=rare 
earth, T=Fe, Ru, Os, X=P, As, Sb) show a variety 
of interesting properties, but their physical mechanisms 
are not well understood. Among them, however, Ce- 
skutterudites show a clear systematic trend that X=P 
and As compounds are insulators whereas X=Sb yield 
metals. 1 This is due to the larger atomic size from P to 
Sb, and hence larger hybridization with f electron on Ce 
in the X12 cage. Similar but weaker trend is seen also for 
T=Fe to Os. In fact, all X=P compounds are insulators 
with a large energy gaps of the order of 0.1 to 0.03 eV, 
whereas CeRu4Sbi2 is a metal at low temperatures ac- 
cording to the transport measurement. 2,3 Applied pres- 
sure renders CeRu4Sbi2 to an insulator for P > 6GPa, 
supporting the above picture. 4 

It is strange that nevertheless CeOs4Sbi2 shows an 
insulating behavior at low temperatures, although the 
band calculation shows that it is a semimetal. 5 Both Ru 
and Os compounds show similar pseudo-gap structures 
and mid-infrared peak in the optical conductivity spec- 
tra 6 ' 7 similar to the Kondo insulators, 8 ' 9 but CeRu4Sbi2 
shows a metallic behavior below 80 K. Therefore, It is an 
interesting issue whethere Ce-skutterudites with insulat- 
ing behaviors at low temperatures are classified into the 
Kondo insulators or not, and what is a mechanism for 
the gap. It was recently found that the insulating be- 
havior of CeOs4Sbi2 may be due to a formation of some 
ordered state (possibly spin density wave (SDW)), 1 to 
which we have proposed a theoretical model and calcula- 
tion explaining the occurence of SDW and the anomalous 
phase diagram in this compounds. 10 Therefore, we focus 
on CeRu4Sbi2 in the present paper. 

We previously proposed a tigh-binding band model 
with full f-electron characters on Ce with the spin-orbit 
(s.o.) interaction and the crystal electric field (CEF), and 
a simplified A u top-most valence band composed of the p 
states on X12 clusters. 11 The calculated optical conduc- 
tivity successfully explained the temperature dependence 



of the mid-infrared (IR) peak and the filling-up of the 
gap with raising temperature by inclusion of the strong 
correlation effects. However, we have assumed a finite 
band gap in the previous work to focus on these temper- 
ature dependence of the (pseudo-)gap structures. Also 
the band calculation gave a finite gap. 5 It is now clear 
that CeRu4Sbi2 is a metal at the ground state, perhaps a 
semimetal, 3 so that we improve our model with a slight 
simplification and calculate the transport properties in 
the present paper. 

In CeRu4Sbi2, the resistivity shows a large peak at 80 
K after subtraction of the data of the La-counterpart. 
Seebeck coefficient is positive and also shows a large 
peak at nearly the same temperature. 2 ' 3 ' 12 Non-Fermi- 
liquid(NFL)-like behaviors are found at very low temper- 
atures below several Kelvin, 2, 12 but we do not take them 
into our scope here. 

In order to explain all these features except NFL, we 
propose a simplified tight-binding model, which captures 
the essential features of the band calculation. Using this 
model and introducing the correlation effect within the 
framework of the dynamical mean field approximation 
and the iterative perturbation theory, 15 the tempera- 
ture dependences of the optical conductivity, resistiv- 
ity and the Seebeck coefficient are calculated, which can 
explain the experiments rather well. Theoretical under- 
standing of the transport properties of the skutterudite 
compounds is very important since they have attracted 
much interest because of its high potentiality for an ap- 
plication as a new good thermoelectric device. 13 ' 14 

The basic crystal structure of filled skutterudites is 
the body-centered cubic (bec) lattice composed of the 
X12 clusters filled by R in the center. The f orbitals of 
rare-earth atoms strongly hybridize with the p orbitals 
on X12 clusters close the Fermi level. The X12 cluster has 
the Th symmetry, but its effect on J = 5/2 f 1 state in Ce 
is the same as that of Oh- Therefore, we use the notation 
for Oh in the fallowings. 
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The present band model consists of the single wide 
band (called c-band hereafter): 



k z a 



e k = t\ cos —j^ - cos cos 

+t 2 (cos k x a + cos k y a + cos k z a) (1) 

with the nearest and the next-nearest neighbor hoppings 
on the bec lattice, which represents the top-most p band 
of X12 clusters. The lattice constenat a is set equal to 
unity in the following. For f-electrons, the neutron scat- 
tering experiment cannot resolve the CEF for the mo- 
ment. 16 We here assume that the the ground state of f 
electron is T 7 doublet as in CeOs 4 Sb 12 , 17 and take the 
doubly degenerate f-band, but since CEF is not well sep- 
arated, it might be the case that the 6-fold J=5/2 states 
could be a better starting model. We neglect this prob- 
lem and take the simplest model, but the degeneracy may 
change the following results only quantitatively. Then, 
we assume the following simple f-band: e£ = Ef + as^, 
where a is a small factor representing the width of the f- 
band dispersion. This would be a reasonable choice since 
both X 12 and R occupy the same bec sites. These two c 
and f bands hybridize with each other. 
Then Hamiltonian is given by 



H 
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e k C k<j C k<r 



E £ k/k<j-/' ko ' 
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E 
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where cj c<r (/^ (T ) creates a c (f) electron with the mo- 
mentum k and spin a (representing T? doublet), and 
f t 

n\ a =j ia fia with the site index i. V is the hybridization 
between conduction and f electrons. Its k-dependence is 
neglected. These choices in our model band immensely 
facilitate the calculations of the strong correlation ef- 
fects, which will be seen later. Actually, of course, the 
f-dispersion is created through the X12 clusters, so that 
it must be self-consistently determined with e£ and the 
hybridization V^. U is the Coulomb repulsion between f 
electrons. 

Without the Coulomb interaction, the diagonalized en- 
ergy dispersions become = [e£ + e£ ± {(e£ — e£) 2 + 
W 2 }^ 2 ]/2, which are the so-called hybridized bands, 

£ 

but since e£ has a finite dispersion, the two bands can 
have finite overlap if a is not too small. 

Considering the empirical trend of the lattice con- 
stant and the energy gap of Ce-skutterudites, CeRu4Sbi 2 
must be a semimetal. 1 Our assumed band structure with 
ti=1.0, t 2 = 0.1, a=0.02, E { =0.3 and V=0.1 is shown 
in Fig. 1, which has a semi- metallic character: the top of 
lower band (r point) and the bottom of upper band (H 
point) slightly touch the Fermi level, respectively, and 
overlap by about 0.03. Although our band structure is 
quite simplified, Fig. 1 captures the low-energy structure 
of the band calculation for CcRu4Sbi2 18 rather well. The 
density of states (DOS) p c (s) for e£ and those for E^ are 
shown in Fig. 2. 

For treating the correlation effects, we use the dynam- 
ical mean-field theory and the iterative perturbation the- 




P N H 



Fig. 1. The model band structure of CcRu4Sbi2 calculated with 
the parameters ti=1.0, t 2 = 0.1, a=0.02, E t =0.3 and V=Q.l is 
shown. The dotted line represents the Fermi level. 




Fig. 2. The density of states p c (e) for e£ (solid line) and those 
for (dotted lines) are shown. 



ory, 15 which is a second-order perturbation theory mod- 
ified so as to interpolate between the weak and strong 
correlation limits. Because of our choice e£ = Ef + ae^. 
and the fc-independent V, the fc-summation is converted 
into a single energy integration on e£ with the density of 
states for this band, p c (e). Therefore, the fc-summation 
over the Brillouine zone is necessary only once at the 
begining and for the calculation of the dynamical con- 
ductivity given by 20 



Rea(u 



2^e 2 C^E 



de 



f(e)-f(e + u) 



Pk(e)pk(e + w) + ap c 1 /{e)p c J{s + u) 



+a 2 p[(e)pi(e + w] 



(3) 



Here we neglected the vertex correction and replace the 
velocity matrix elements by an average v 2 cx except its 
scaling factor a, /(e) the Fermi function and p^(e) = 
(— l/7r)ImG^(e + i5) (7 =c or f or cf) with the Green 
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functions 

Gt{e) = 

G c J(s) = 

G f k (e) = 



e - el - 



V 2 
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When the imaginary parts of the self-energy and a are 
small, Reer(aj), are dominated by the current carried by 
the c-electrons and can be approximated by 21 



Reer(a 



2ire 2 v 2 ]T 



xlm- 



W + E*( e £,fc)-E^(4,fe) : 
which is not a bad approximation. Here, 

V 2 



£ c (e,fc) 



(7) 



(8) 



and R (A) denotes the retarded (advanced) function. 

The resistivity p(T) and the Seebeck coefficient S{T) 22 
for correlated systems are given by 



l/p(T) = 2e- 



J deL(e) 



del' 



(9) 



S(T) = --LJ d£ L(e ){ s-p)(- d £ 

I J deL(e) , (10) 

where L(e) = (irN)' 1 J2w( v kx) 2 [I m G£(e, k)] 2 when the 
c-electron dominates and the vertex corrections can be 
neglected. It can be further approximated by L{e) ~ 
p c (e)v 2 x T c (e) similar to eq.(7) with the relaxation time 
of c electrons: Tc(e^)^ 1 ~ — 2ImE c (e k , k). Both approxi- 
mations yield the similar results, as will be shown below. 



Fig. 4. The experimental 6 and theoretical results for the dynam- 
ical conductivity for various temperatures are shown. 



We numerically solve the f-electron self-energy E/(e) 
for [7=0.3 by iteration, and calculate the above- 
mentioned quantities. The fc-summations are carried over 
400 3 mesh points in the 1/8 of the extended cubic Bril- 
louin zone. Because of the correlation effects, the quasi- 
particle DOS becomes temperature-dependent, and the 
sharp Kondo peak grows up on the Fermi level at low 
temperatures in addition to the upper Hubbard peak 
(Fig. 3). The appearance of the Kondo peak enhances the 
temperature-dependence of the physical quantities com- 
pared to the rigid band calculation. 

Choosing the scale as t\ = 7000 K, the calculated op- 
tical conductivity (by the correct formula, eq.(3), but 
with only the cc-part) in Fig. 4 shows a mid-IR peak at 
0.1 eV in accord with the experiment, but the width 
is narrower. There is a large pseudogap below the mid- 
IR peak which is due to the direct (fc-conserving) inter- 
band transition, but a Drude peak at low temperature 
and low frequencies is seen since it is a semimetal. The 
main descripancy is that the experimental spectrum at 
the lowest temperature has a long tail below the mid- 
IR peak, which may be due to the indirect transitions, 
as was pointed out recently by Okamura, et al. and the 
present authorfor the Kondo insulator YbBi2. 8,9 ' n The 
temperature-dependence is weaker than the experiments, 
but the psuedo-gap of the order of 0.1 eV is filled up 
already at T = 300K, which is due to the many-body 
effect. 9,11 Despite some descripancies, the overall struc- 
tures and the tepmerature-dependence are rather well 
reproduced by the present simple model. 

The resistivity and the Seebeck coefficient with and 
without the relaxation approximation are compared with 
the experiment in Figs. 5 and 6. The calculated peak po- 
sition agrees with the experiment in the Seebeck coeffi- 
cient, but is twice higher in the resistivity. Nevertheless, 
our results capture the global features rather well, except 
the strange steep rise of S(T) at very low temperatures 
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f band) to describe the low energy state of CeRu 4 Sbi 2 , 
and calculated the optical conductivity, the resistiv- 
ity and the Seebeck coefficient, taking account of the 
Coulomb interaction between f electrons via the dynam- 
ical mean-field approximation and the iterative pertur- 
bation theory. The results can explain the experimental 
data semi-quantitatively. Improvement of the model and 
the consideration of the non-Fermi-liquid-like behaviors 
found at low temperatures remain to be done in the fu- 
ture. 
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sions. This work was supported by The Grant-in-Aid 
from the Ministry of Education, Science and Culture: 
"Evolution of New Quantum Phenomena Realized in the 
Filled Skutterudite Structure", No. 16037204. 



Fig. 5. The calculated resistivity (the broken line by the relax- 
ation time approximation) is compared with the experimental 
value is shown. (The data of the coressponding La-compound is 
subtracted from the experimental data. 3 ) 



200 




100 200 
T[K] 



300 



Fig. 6. The calculation for the Seebeck coefficient (the broken 
line by the relaxation time approximation) is compared with the 
experiment. 3 



and the shoulder at 20 K. The difference in p(T) at low 
temperatures may be due to the residual resistivity. 

We tried to fit all these data by the present model by 
changing the band parameters and U , but the present 
ones are the best. The remaining discrepancies may be 
due to the oversimplification of the present band model, 
improvement of which is left for a future study. 

In summary, we have constructed the effective tight- 
binding model (periodic Anderson model with dispersive 
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